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Abstract 

We develop and significantly generalize the effective worldvolume theory for higher- 
dimensional black holes recently proposed by the authors. The theory, which regards 
the black hole as a black brane curved into a submanifold of a background spacetime 
— a blackfold — , can be formulated in terms of an effective fluid that lives on a dynamical 
worldvolume. Thus the blackfold equations split into intrinsic (fluid-dynamical) equations, 
and extrinsic (generalized geodesic embedding) equations. The intrinsic equations can be 
easily solved for equilibrium configurations, thus providing an efficient formalism for the 
approximate construction of novel stationary black holes. Furthermore, it is possible 
to study time evolution. In particular, the long-wavelength component of the Gregory- 
Laflamme instability of black branes is obtained as a sound-mode instability of the effective 
fluid. We also discuss action principles, connections to black hole thermodynamics, and 
other consequences and possible extensions of the approach. Finally, we outline how the 
fluid/ AdS-gravity correspondence is related to this formalism. 
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1 Introduction 



In recent work we have identified the origin of the rich variety of higher-dimensional black 
holes in the possibility of having horizons that are much longer along some directions 
than in others pp. In such cases there are at least two different horizon length-scales, and 
we have proposed an effective theory that captures the long-distance physics when these 
scales are widely separate. Focusing on the simplest, universal case of neutral, vacuum 
black holes, the two length scales are associated with the mass and angular momentuirj^], 

In a four-dimensional black hole, by virtue of the Kerr bound J < GM 2 these lengths are 
always parametrically similar, but higher-dimensional black holes (including also black 
rings) with lj 3> £m are known to exist in ultra-spinning regimes where the angular 
momentum for a given mass can be arbitrarily high [21E1E1E]- It appears that essentially 
all their novel features, compared to their four-dimensional cousins, arise from the ability 
to separate these two lengths. This suggests that higher-dimensional black holes must be 
organized according to a hierarchy of scales: 

1- black holes behave qualitatively similarly to the four-dimensional Kerr 

black hole. 

2. I j ~ £m' threshold of new black hole dynamics. 

3. fj » £m- the separation of scales suggests an effective description of long- wavelength 
dynamics. 

The first and second regimes fully involve the non-linearities of General Relativity, but 

in the first case we have no hints of qualitatively new properties of black holes compared 

i 

to four-dimensional ones. In particular we conjecture that for J < J cr j t = aoM(GM) D -s 
(with old a yet-undetermined numerical constant of order one) , the Myers- Perry black holes 
are dynamically stable and unique among solutions with connected regular horizons^]. In 
contrast, as the two scales begin to diverge in the second regime, we have good evidence 
of the onset of new phenomena: horizon instabilities, inhomogeneous ('pinched') phases, 
non-spherical horizon topologies, and absence of uniqueness [H El El [7]. This regime seems 
hard to investigate by means of exact analytical techniques, but on the other hand the 
presence of only one scale in the problem is actually convenient for numerical investigation, 
since one does not require high precision over widely different scales. In the third regime, 
which is the focus of this paper, the existence of a small parameter ImI^j allows the 

1 J = (£\ jf) 1//2 aggregates the effect of all possible angular momenta. 

2 This extends a previous conjecture in [3] about stability of MP black holes to involve their uniqueness 
too. Here J cr it refers to the minimum J for which either instability or non-uniqueness first appear. 
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introduction of efficient approximate analytical methods. In this way, this organization in 
scaled provides an outline for a program to investigate higher-dimensional black holefl 

We base the effective description of black holes in regime 3 on the idea that the limit 
@m/£j of the black hole, keeping its horizon size finite, results in a flat, infinitely 
extended black brane, possibly boosted along some of its worldvolume directions. This 
has been observed in all exact solutions known so far: the horizon of Myers-Perry black 
holes pancakes along the plane of rotation and in the limit becomes a black brane [3], and 
five-dimensional black rings become thin and locally approach the geometry of boosted 
black strings |114 112j . This allows us to identify the variables for the effective description. 
Essentially, the black hole will be regarded as a black brane whose worldvolume spans a 
curved submani/o/rf of a background spacetime — what we refer to as a blackfold. The 
simplest example is the characterization of a (thin) black ring as a circular boosted black 
string, which was worked out in detail in [6] and expanded upon in [131114] . But, as shown 
in PQ, ultraspinning Myers-Perry black holes, as well as a number of other new black holes, 
can also be appropriately captured with this approach. 

Thus, the effective theory of black holes when £m/^j *C 1 is a theory that describes 
how to bend the worldvolume of a black brane in a background spacetime. In this regard 
we treat black branes in a manner similar to other familiar extended objects such as cosmic 
strings or D-branes. The main novelty is that black branes possess black hole horizons, 
so when their worldvolume is spatially compact we obtain a black hole with finite horizon 
area. 

To leading order in the expansion in Im/£j the backreaction of the blackfold on the 
geometry is neglected — it is a 'test' blackfold. Corrections to the geometry are found 
by first computing the linearized gravitational backreaction on the background, which 
is of order (£m/£j) D ~ 3 , then this correction induces a perturbation of the near-horizon 
geometry, and so forth. This can be systematically pursued in the form of a matched 
asymptotic expansion. For higher-dimensional black rings these first corrections have 
been computed [6]. In this paper, however, we remain at the 'test blackfold' level of 
approximation. This is enough to reveal new kinds of black holes, compute their physical 
properties, and also study time-dependent situations and stability. 

Ref. [Ij gave a basic outline of the theory of blackfolds. Here we develop its concep- 
tual basis and considerably improve and generalize the presentation. Its application to 
specific new classes of higher-dimensional black holes will be discussed in a forthcoming 
publication. 

We draw heavily from the beautiful theory of classical brane dynamics developed by 
Carter in [15^1 . In this respect, we may say that a main part of our contribution is to 

3 The approach to extremality in a black hole introduces a long length scale transverse to the horizon 
and allows to decouple a different sector of the physics, namely the near-horizon region [8]. 

4 See [9] QJj] for brief reviews of higher-dimensional black holes and [5] for a more extensive one. 
5 Our title deliberately highlights this. 
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apply this theory to black branes, and to interpret the results in the context of higher- 
dimensional black hole physics according to the general considerations above. But we 
also emphasize that the general classical brane dynamics, regarded as a long- wavelength 
effective theory, must take the form of the dynamics of a fluid that lives on a dynamical 
worldvolume. Black branes correspond to a specific type of fluid, with a certain equation of 
state and with specific values for transport coefficients. At the leading order that we work 
in this paper, the fluid is a perfect one and the brane equations are the Euler equations 
for the fluid — intrinsic equations — plus a generalization of the geodesic equation for the 
motion of a p-brane — extrinsic equations for the worldvolume embedding. We believe 
that in principle it should be possible to incorporate higher-derivative corrections to these 
and compute transport coefficients by performing a derivative expansion of the underlying 
microscopic theory, in this case Einstein's theory. 

Closely related precedents of a mapping of black hole dynamics to fluid dynamics are 
the 'membrane paradigm' |16j . and the more recent 'fluid/AdS-gravity correspondence' 
[T7] . As we shall argue near the end, the fluid/ AdS-gravity correspondence can be embed- 
ded within the approach we advocate here — in fact it has been an important influence 
in developing it. The general arguments discussed above indicate that a fluid-dynamical 
description should indeed be expected to exist for any long-wavelength fluctuations around 
an equilibrium state. From this perspective, perhaps the main qualitative novelty of our 
approach is that the existence of a hierarchy of scales in higher-dimensional black holes 
makes this effective theory useful not only for studying fluctuations, but also for construct- 
ing and analyzing in a very general manner novel kinds of stationary (i.e., equilibrium) 
black holes, including vacuum solutions. 

The outline of the paper is the following: Section [2] develops the conceptual basis 
underlying the blackfold approach as a worldvolume theory of the dynamics of black 
branes. Section [3] presents a main result of this paper: the blackfold equations, a set of 
coupled non-linear differential equations for the collective coordinates of a neutral black 
brane. Section [J] focuses on the important case of stationary blackfolds, for which the 
intrinsic subset of these equations can be explicitly solved. Section [5] analyzes the issues 
raised by the possible presence of boundaries of the blackfold worldvolume. In section [6] 
we describe how to compute the physical magnitudes of a blackfold. Section [7] presents 
an action principle for stationary blackfolds. This is useful for practical calculations, but 
also admits a simple and appealing interpretation in terms of black hole thermodynamics. 
Section [8] discusses briefly the stability of blackfolds exhibiting how the approach can 
uncover in a remarkably simple way the Gregory-Laflamme instability of black branes. 
We close in Section [9] with a discussion of the relation of blackfolds to other effective 
theories of black hole dynamics, in particular the fluid/AdS-gravity correspondence. In the 
appendix we collect a number of technical results on the extrinsic geometry of submanifold 
embeddings. 
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Notation and terminology: 

For clarity and later reference, we summarize here some of our notation. 

For a blackfold of p spatial dimensions in D-dimensional spacetime it is convenient to 
introduce 



The codimension of the blackfold worldvolume is n + 2. 

Spacetime (background) and worldvolume magnitudes are denoted and distinguished 
as follows: 

• Spacetime coordinates (and embedding functions): X 11 , fi, v . . . = 0, . . . , D — 1. 
Background metric: g^ u . 

Background metric connection: r^„. 
Background covariant derivative: V^. 

• Worldvolume coordinates: a a , a, b . . . = 0, . . . , p. 
(Induced) worldvolume metric: 7^. 
Worldvolume metric connection: { 6 ° }• 
Worldvolume covariant derivative: D a . 

Indices /x, u, . . . are lowered and raised with g^, indices a,b, . . . with j a j,. 

We use the same letter for a background tensor tangent to the worldvolume, t^... 1 '"', and 
for its pullback onto the worldvolume, t a ,, b "' (the only exception is the first fundamental 
form h^y and the induced metric j a b). 

f2( n ) denotes the volume of the unit n-sphere. Jlj denotes the angular velocity of the 
blackfold in the i-th. direction. 

V(p) is the volume of a spatial section of the blackfold. Vi is the spatial velocity field 



Note that we refer to 'long- wavelength' and not 'low-energy', effective theory, the 
reason being that in classical gravity, without h, such notions are not equivalent (indeed 
large energies typically imply long distances in classical gravity). 

2 Effective worldvolume theory 

We present the effective theory of blackfolds trying to highlight the similarities with the 
field-theoretical effective description of other extended objects, such as cosmic strings or 
D-branes. The main differences with these are, first, that the short-distance degrees of 
freedom that are integrated out are not those of an Abelian Higgs model nor massive 
string modes, but rather purely gravitational degrees of freedom. Second, the extended 
objects — curved black branes — possess black hole horizons. We obtain the equations 
using general symmetry and conservation considerations, rather than doing a detailed 
derivation from first principles. 



n = D — p — 3 . 



(1.2) 



on the worldvolume, and V = v/^j ■ 
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2.1 Collective coordinates for a black brane 

Schematically, the degrees of freedom of General Relativity are split into long and short 
wavelength components, 

^ = {^° ns) >#° rt) }- (2-1) 



The Einstein-Hilbert action is then approximated ag, 

lEH= \tcS d ° x ^R ~ J^Q J d D x^-g(^)R(^ S ) + I es[g ^)^] , (2.2) 

where Iesio^u^, 4 1 ] 1S an effective action obtained after integrating-out the short-wavelength 
gravitational degrees of freedom (precisely what we mean by this will be made clear in 
sec. I2.2p . The coupling of these to the long- wavelength component of the gravitational 
field is captured through a set of 'collective coordinates' that we denote schematically by 
4>. Our first task is to identify these effective field variables and the length scales that 
allow this splitting of degrees of freedom. 

The main clue to the nature of the effective theory comes from the observation that the 
limit O-m/O-j — * of known black holes, when it exists, results in flat black branes. Thus 
we shall take the effective theory to describe the collective dynamics of a black p-brane. 
Its geometry inD = 3+ p + n spacetime dimensions is 



^°p- brane ' 



- ) dt* + J>^) 2 + + r*dnl +1 . (2.3) 

' i=l 1 r n 



The coordinates a a = (t,z l ) span the brane worldvolume. A more general form of the 
metric is obtained by boosting it along the worldvolume. If the velocity field is u a , with 
u a u b r]ab = — 1 then 

Vab + ^u a u b J da a da b + ^ + r 2 dn 2 n+1 . (2.4) 

The parameters of this black brane solution consist of the 'horizon thickness' Tq, the 
p independent components of the velocity u (say, its spatial components u l ), and the 
D — p — 1 coordinates that parametrize the position of the brane in directions transverse 
to the worldvolume, which we denote collectively by X-^. The D collective coordinates of 
the black brane are 

<f>(a a ) = {X ± (a a ),r (a a ),u i (a a )} (2.5) 

and in the long- wavelength effective theory one allows dX- 1 -, lnro and U{ to vary slowly 
along the worldvolume, W p +i, over a length scale R much longer than the size-scale of the 
black brane, 

R > r . (2.6) 



6 For clarity of presentation, at this initial stage we consider that both long and short degrees of freedom 
obey vacuum gravity dynamics, i? M „ = 0, but this can be easily generalized, see sec. 13.51 below. 
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Typically the scale R is set by the smallest intrinsic or extrinsic curvature radius of the 
world volume. Observe that we require slow variations of dX^~, not of X-^. Like the 
longitudinal velocities u a , the transverse 'velocities' dX 1 - can be arbitrary. 

In order to preserve manifest diffeomorphism invariance it is convenient to introduce 
some gauge redundancy and enlarge the set of embedding coordinates of the worldvolume 
of the black brane to include all the spacetime coordinates X^(a a ). From this embedding 
we can compute an induced metric 

la b = g^ g) d a X»d b X» . (2.7) 

This is naturally interpreted as the geometry induced on the worldvolume of the brane. 
To understand what this means, regard the split between degrees of freedom as follows: 
the long- wavelength degrees of freedom live in a 'far-zone' r > ro, and they describe the 
background geometry in which the (thin) brane lives. Then (|2.7p is the metric induced on 
the brane worldvolume. The short-wavelength degrees of freedom live in the 'near-zone' 
r <C R. In the strict limit where R —* oo, the near- zone solution is (|2.4|) . but when R is 
large but finite, the collective coordinates depend on a. Also, the long and short degrees 
of freedom interact together in the 'overlap' or 'matching-zone' ro <C r <C R, where the 
metrics g^° ng ^ and g^v° must match. Then the near-zone metric for the black brane 
must be of the form 

/ r n (rr) \ dr 2 

Hbort) = + ^T^)u h {a) J da a da h + + r 2 d^ +1 + . . . . ( 2 . 8 ) 

The dots here indicate that, without additional terms, in general this is not a solution 
to the Einstein equations. These equations contain terms with gradients of lnro, u a and 
7 a 6. However these terms can be seen to come multiplied by powers of ro so they are 
small when tq/R -C 1. Then we can consider an expansion of the equations in derivatives 
and add a correction to (12. 8p to find a solution to the Einstein equations to first order in 
the derivative expansion. A subset of the resulting Einstein equations can be rewritten 
as equations on the collective field variables cp(cr). An important requirement is that the 
perturbations preserve the regularity of the horizon, and to this effect working in a set 
of coordinates (Eddington-Finkelstein type) different than the ones above may be more 
appropriate. 

The development of this line of argument, which can be regarded as a blend of the ideas 
for the effective descriptions of black hole dynamics in |18|, [19] (and references therein), 
and in [T7], produces a systematic derivation of the blackfold equations. This is however 
a technically involved approach that we hope to discuss elsewhere. Here we shall instead 
follow a less rigorous but quicker and physically well-motivated path, relying on general 
effective-theory-type of arguments that allow us to readily obtain the blackfold formalism 
valid to lowest order in the derivative expansion. As we will see, this is the 'perfect fluid' 
and 'generalized geodesic' approximation. The more systematic method outlined above 
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would be needed to go beyond these approximations and account for dissipation and effects 
of internal structure and gravitational self-forc^E 

We have considered a black brane in (|2,3p , (|2.4p that is not rotating along the transverse 
(n + l)-sphere, nor have we included any possible deformations of it. This is just a 
simplification and is not essential. Note first that ultraspins in this sphere, with rotation 
parameter a oc J/M 3> ro, can indeed be considered in the blackfold approach, but then 
the starting point must be the black brane limit that results when a/ro — > oo. On the 
other hand small spins, with a < ro, can also be included easily to the order that we work 
in this paper. The reason is that the modifications to the blackfold dynamics introduced 
by these spins only enter at a higher order in the expansion in ro/R. This is familiar, for 
instance, in that spin effects on the worldline of a test particle enter through couplings to 
the background curvature tensor: they reflect the internal structure of the particle, which 
clearly is a higher-order correction. Thus the dynamics associated to internal spin and 
polarization effects of the spheres of size ro is effectively integrated-over without affecting 
the lowest-order formalism. It must be noted, though, that higher-dimensional neutral 
black holes exhibit zero-mode deformations at discrete values of the spin when n > 3 
(when n = 1,2 such deformations appear to be always massive modes) |7|. For these 
spins, these deformations can be excited at arbitrarily low frequencies and therefore must 
be added to the collective coordinates of the black brane. Other than this, the internal 
spin can be treated just like a conserved charge in the worldvolume. As such, it must 
satisfy a continuity equation but, since the gravitational effects of this spin fall off faster 
at large r than those of the mass, the computation of the effective worldvolume stress 
tensor in the next subsection is unaffected. The blackfold equations also apply in the form 
below in the presence of internal spins. 

2.2 Effective stress tensor 

By the phrase 'integrating out the short-distance dynamics' we mean that the Einstein 
equations are solved at distances r <C R and then the effects of the solution at distances 
r 3> ro are encoded in a stress-energy tensor that depends only on the collective coor- 
dinates. The stress tensor is such that its effect on the long-wavelength field gjt° ng ^ is 
the same as that of the black brane at distances r > ro. For reasons that will become 
apparent as we proceed, it is both simpler and more convenient to work with an effective 
stress-energy tensor rather than with an effective action. In any case, nothing is lost since 
we work at the classical level. 

The effective equations from (|2.2p are 

<° ng) " ^ (1 ° ng) rfr g) = 8nGT^ , (2.9) 

7 Also, see [20] for an approach to the calculation of these corrections (for zero-branes) in an effective- 
theory framework somewhat akin to the spirit in this paper. 
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where the effective worldvolume stress tensor is 



T cfT 2 ^cff 
I X ^ v 



(2.10) 

W P+ i 



We now argue that the appropriate notion for this effective stress-tensor that captures 
the coupling of the short- wavelength degrees of freedom to the long- wavelength ones, is 
the quasilocal stress-energy tensor introduced by Brown and York [21) . This is defined by 
considering a timelike hypersurface that lies away from the black brane and encloses it 
by extending along the worldvolume directions and the angular directions £lr n +u, i.e., the 
hypersurface acts as a boundary. Actually, as we explained above, the angular directions 
are integrated over in our description (and to leading order they do not play any role), 
so we can simplify the discussion by focusing exclusively on the worldvolume directions 
of the boundary. If the boundary metric (along worldvolume directions) is j a b then the 
quasilocal stress tensor is 

rp (quasilocal) _ 2 SI d 

1 ab ~ .^R^ab ' y l - ll > 



-7 5^° 

where I c i is the classical on-shell action of the solution. For our purposes, this is the action 
where the short-distance gravitational degrees of freedom, r <C R, are integrated and so 
it must be the same function of the collective variables as I e ft- Together with the relation 
fl22D this implies that we can identify §LWj) with (|2TTT|) . 

It is shown in [21] that the Einstein equations with an index orthogonal to the boundary 
are first-order equations equivalent to the equation of conservation of the quasilocal stress 
tensor, 

^(quasilocal) =0, (2-12) 

where D a is the covariant derivative associated to the boundary metric 7 a ft. Hence, solving 
the equations (|2.12p is equivalent to solving (a subset of) the Einstein equations. 

Since we identify the stress tensors (|2.10p and (|2.1ip . henceforth we drop the super- 
scripts from them. We also drop the superscript ( lon s) from the background metric g^ v . 

The effective stress tensor is computed in the zone ro <C r <C R, where the gravitational 
field is weak and the quasilocal stress tensor T ab is, to leading order in tq/R, the same as 
the ADM stress tensor. For the boosted black p-brane (|2.4p one can readily compute it 
and find 



T ab = (nu a u b - V ab ) . (2.13) 

After introducing a slow variation of the collective coordinates the stress tensor becomes 

T ab (a) = ^k(cr) (nu a (e)u\e) - ^(a)) + . . . (2.14) 

where the dots stand for terms with gradients of lnro, u a , and 7 a fe, which we are taking 
to be small and are neglected in this paper. 
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2.3 General branes: fluid perspective 

On general grounds, the long- wavelength effective theory for any kind of brane will take 
the form of a derivative expansion for an effective stress-energy tensor that satisfies the 
conservation equations (|2.12p . This is the dynamics of an effective fluid that lives on the 
worldvolume spanned by the brane. If the worldvolume theory is isotropic, then to lowest 
derivative order the stress tensor is that of an isotropic perfect fluid, 

T ab = (e + P)u a u b + Pj ab , (2.15) 

with energy density e, pressure P and velocity u a satisfying 

u a u b lab = -l. (2.16) 

Thermodynamics provides the universal macroscopic description of equilibrium configura- 
tions, and fluid dynamics is the general long- wavelength description of fluctuations under 
the assumption of local equilibrium. So in general there will be an equation of state, which 
we write in the form P(e), and the system will obey locally the laws of thermodynamics 

de = Tds (2.17) 

and Euler-Gibbs-Duhem relation 

e + P = Ts (2.18) 

where T is the local temperature and s the entropy density of the fluid in its rest frame. 
The fluid may also carry additional conserved charges, but we do not consider these in 
this paper. 

For a black brane, (|2.14p tells us that the effective fluid has 

6 = %±£(f> + l)rff, P = —L_ £ . (2.19) 
167rG n+1 

Moreover, in the rest frame of the fluid the Bekenstein-Hawking identification between 
horizon area and entropy 

- = ^V- 1 ( 2 - 2 °) 
and between surface gravity and temperature 

r = - — 2.2i 

4irr 

is well known to reproduce the correct thermodynamic relations (|2.17|) . f|2. 18(1 . 

Going beyond the perfect fluid approximation (|2.15p . the stress tensor will acquire 
dissipative terms proportional to gradients of lnro, u a , j ab . As discussed above, these are 
neglected in this paper. At any rate their effects are absent for stationary configurations. 



10 



3 Blackfold dynamics 



We have argued that the general effective theory of classical brane dynamics can be formu- 
lated as a theory of a fluid on a dynamical worldvolume. The fluid variables must satisfy 
the intrinsic equations (|2.12p . and they will be coupled to the 'extrinsic' equations for the 
dynamics of the worldvolume geometry, which we still have to determine. To this effect, 
in the next subsection we introduce a few notions about the geometry of worldvolume 
embeddings. More details and proofs are provided in the appendix. 

3.1 Embedding and worldvolume geometry 

Given the induced metric on W p +i, (|2.7|) . the first fundamental form of the submanifold 
is 

hr = d a x^d b x u j ab . (3.1) 

Indices fj, , v are raised and lowered with g^ u , and a, b with j ab . Defining 

-L^ti/ = Qfj,v hpi/ (3-2) 

it is easy to see that the tensor h^ u acts as a projector onto W p +i, and J-^ u along directions 
orthogonal to Wp+i. 

Background tensors with support on W p +i can be converted into worldvolume 

tensors t a "5... and viceversa using d a X^. For instance, the velocity field 

vP = d a X» u a , (3.3) 

preserves its negative-unit norm under this mapping. 

The covariant differentiation of tensors that live in the worldvolume is well defined 
only along tangential directions, which we denote by an overbar, 

V M = y V„ • (3.4) 

Note that in general Vpt^"V... h as both orthogonal and tangential components. The 
tangentially projected part is essentially the same as the worldvolume covariant derivative 
D c t a "b... f° r the metric 7 a b, both tensors being related via the pull-back map d a X^. In 
particular, the divergence of the stress-energy tensor 

T» v = d a X^d b X u T ab (3.5) 

satisfies (see (|A.2ip ) 

h p u V^ u = d b X p D a T ab . (3.6) 
The extrinsic curvature tensor 

K^ p = h^V v K p (3.7) 
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is tangent to W p +i along its (symmetric) lower indices /jl, u, and orthogonal to W p+ i along 
p. Its trace is the mean curvature vector 

k p = h^K^o = V^ p . (3.8) 

Explicit expressions for the extrinsic curvature tensor in terms of the embedding functions 
X^(a a ) can be found in the appendix. 

3.2 Blackfold equations 

The general extrinsic dynamics of a brane has been analyzed by Carter in [15]. The 
equations are formulated in terms of a stress-energy tensor with support on the p + 1- 
dimensional worldvolume W p +i satisfying the tangentiality condition 

j_p^ = o . (3.9) 

The basic assumptions are that (i) this effective stress-energy tensor derives from an 
underlying conservative dynamics (in our case, General Relativity), even if the macroscopic 
(=long- wavelength) dynamics may be dissipative; and that (ii) spacetime diffeomorphism 
invariance holds, or equivalently, the worldvolume theory can be consistently coupled to 
the long- wavelength gravitational field g^v- Under these assumptions, the stress tensor 
must obey the conservation equations 

V M TW = 0. (3.10) 

These are in fact the generic equations of motion for the entire set of worldvolume field 
variables 4>(a a ), both intrinsic and extrinsic: we can decompose (|3. 10[) along directions 
parallel and orthogonal to W p+ i as 

= T^K lxv P + d h XPD a T ah (3.11) 

where in the last line we used (|3.6j) and (|3.7p . Thus the D equations (|3,10p separate into 
D — p—1 equations in directions orthogonal to W p +i and p + l equations parallel to W p +i, 

T^K^yP = (extrinsic equations) , (3-12) 
D a T ab = (intrinsic equations) . (3.13) 

Let us now apply the equations (|3.1U|) onto the generic stress-energy tensor of a perfect 
fluid on the worldvolume 

r = (e + P)uV+PF. (3.14) 

We find 

uWVve + (e + P)( U p + U PV u u u ) + (h^ + u Pu u )V v P + PRP = , (3.15) 
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where 

u = u u V u u (3.16) 

is the acceleration of vf 1 . These are the general equations, to leading order in the derivative 
expansion, for the dynamics of a classical brane with worldvolume spatial isotropy, possibly 
supplemented by conservation equations for charges, if present. A familiar example is a 
Nambu-Goto-Dirac brane, with = — \P\ha V and VP = Ve = Vit = 0, for which the 
extrinsic equations, K p = 0, require that the worldvolume be a minimal submanifold. But 
any classical brane will satisfy equations of this form. 

As usual the projection of (|3.15[) onto u is the continuity equation for the energy of 
the fluid, 

/V„e+(£ + P)V/ = 0, (3.17) 
while the projections orthogonal to u 

(e + P)u» = -(U LU + u p u u )V p P - PK» (3.18) 

say that the force that accelerates an element of the fluid is given along worldvolume di- 
rections by pressure gradients (Euler equation) and in directions transverse to the world- 
volume by the extrinsic curvature. 

For the specific stress tensor of a neutral black brane, ()2. 19[) . the equations f)3. 15|) 
become, after a little manipulation, 

+ — !— ^V,/!/ = -i^ + V M lnr . (3.19) 
n + 1 n 

These blackfold equations describe the general collective dynamics of a neutral black brane. 

Again, we can decompose them into different projections. In directions orthogonal to 
the worldvolume we have 

K p = n± p f,ii p . (3.20) 

The equivalence of this equation to (|3.12p follows by using (|A. 12|) . 
The equations parallel to the worldvolume are 

ha V u v H l —UuSvU v = V M In r , (3.21) 

n + 1 

which we can also write using worldvolume indices and derivatives, 

u a H l —u a D b u b = d a In r , (3.22) 

n + 1 

with ii b = u c D c u b . Thus the temporal and spatial worldvolume gradients of tq determine 
the worldvolume acceleration and expansion of u, respectively. 

Although we have emphasized the fluid-dynamical interpretation of the equations, it is 
interesting to observe that the extrinsic equations (|3.12p . when written explicitly in terms 
of the embedding X p (a a ) become 

T ah u p {d a d b x a + r° u d a x»d b x u ) = o , (3.23) 
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or alternatively 

T ab {D a d b X? + Y^daX^X") = (3.24) 

(see eqs. (|A.28|) and (|A.30P ), These can be regarded as generalizations to p-branes of the 
geodesic equation for free particles, or more simply, of "mass x acceleration= 0". 

Blackfolds differ from other branes in that they represent objects with black hole hori- 
zons. In the long-distance effective theory we lose sight of the horizon, since its thickness 
is of the order of the scale ro that we integrate out. But the presence of the horizon is 
reflected in the effective theory in the existence of an entropy and in the local thermo- 
dynamic equilibrium of the effective fluid. Indeed, we shall assume that the regularity 
of the event horizon under long-wavelength perturbations — including those that bend 
the worldvolume away from the flat geometry or that excite the effective fluid away from 
equilibrium — is satisfied when the blackfold equations, which incorporate in particular 
local thermodynamic equilibrium, are satisfied. A proof of this statement requires the rig- 
orous derivation, outlined at the end of sec. 12.11 of the blackfold equations in the derivative 
expansion of Einstein's equations, which is outside the scope of this paper. However, there 
is already significant evidence that horizon regularity is preserved for blackfold solutions. 
First, analyses of the perturbations of black strings that bend them into a circle [6] (and 
extensions thereof to branes curved into tori) show that the extrinsic equations (13.121) 
are equivalent to demanding absence of singularities on or outside the horizon. Second, 
the intrinsic, hydro dynamical perturbations of a black brane in AdS have been studied 
in detail in [17] and shown to be consistent with horizon regularity. Note however that 
the regular solution to higher orders may not preserve the same symmetries as the lowest- 
order solution. In particular, in some cases horizon regularity may require to abandon 
stationarity of the configuration at higher orders. In most instances it appears easy to de- 
cide from the physics of the problem whether such an effect is expected (specific examples 
will be discussed elsewhere), but it would be good if precise conditions could be stated in 
generality^ 

3.3 The metric at all length scales: Matched asymptotic expansion 

Under the splitting in (|2.2p . the set of field variables in the system are the collective 
worldvolume fields, intrinsic and extrinsic, and the background gravitational field g^ v . 
The complete set of equations are the extrinsic equations f|3. 12|) , intrinsic equations f|3.13j) , 
and backreaction equations (|2.9|) . Since they are a consequence of general symmetry and 
conservation principles, these equations retain their form at any perturbative orderjj The 

8 The example studied in [22] appears to fall outside the remit of our approach since the blackfold (a 
small black hole) is not a small perturbation of the background spacetime. 

9 This, however, is a somewhat formal statement due to the appearance of gravitational self-force di- 
vergences on the worldvolume that must be dealt with carefully [18| . Ref. |23j shows how the equation of 
stress tensor conservation can be used as the basis to obtain these corrections to particle motion. 
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specific form of the stress tensor, as well as the background metric, will in general be 
corrected at higher orders. 

The only equations that one has to solve at the lowest order are those that suffice 
to ensure that can be consistently coupled to the long-wavelength gravitational field. 
These are just the intrinsic and extrinsic equations, and backreaction is neglected. The 
explicit blackfold equations (|3.19p that result are valid only for test branes. 

From the point of view of effective field theory one is interested only in quantities 
that are measured in the long-wavelength regime. The short-wavelength dynamics enters 
only to determine the coefficients in the effective stress tensor, which can be computed by, 
e.g., matching the calculations of some observables [20]. However, in General Relativity 
one is often interested in also having an explicit solution, even if an approximate one, for 
the geometry at all scales, including near the horizon of the black hole. The systematic 
way to construct an explicit metric in an expansion in ro/R is through the method of 
matched asymptotic expansions [18J. In the context of blackfolds this was discussed in [6] 
(following |24| I25j). and explicitly applied to the construction of higher-dimensional black 
rings. We review it briefly here. 

As described in sec. 12.1} the full geometry splits into near- and far-zones, that share a 
common overlap-zone. To zeroth order in r$/R the near-zone metric is (|2.4p and the far- 
zone metric is the background metric g^ u . This is as far as we go in this paper in terms of 
providing explicit solutions to the Einstein equations: this is the test-brane approximation. 
The next order involves the gravitational backreaction of the brane: the equations (|2.9|) 
are linearized around the background and solved with the distributional worldvolume 
source T^ u . The blackfold equations would appear here among the Einstein equations 
as constraints (first-order equations), but we can assume they have already been solved 
at zeroth order. The solution of the linearized equations with appropriate asymptotics 
(e.g., asymptotic flatness) produces a corrected far-zone metric with corrections of order 
(ro/R) n . Its value in the overlap-zone provides new asymptotic conditions for the near- 
zone solution. The next step is to perturb the metric (|2,4p linearly, with the boundary 
conditions that the horizon remains regular and that in the overlap zone the metric matches 
the corrected far-zone solution. In this manner we produce a new, corrected solution, at 
all scales. This process can then be iterated to higher orders. 

3.4 Lumpy blackfolds 

The formalism developed in this paper can easily be applied to other branes once their 
effective equation of state is known. In particular, there are other neutral black branes 
in vacuum gravity than the 'smooth' black branes of f|2.3|) : refs. [26\ I27j have shown that 
'lumpy' black branes exist, branching off at the threshold of the Gregory-Laflamme (GL) 
instability [28]. Their horizons are inhomogeneous on a scale ~ ro, so this small-scale 
inhomogeneity is averaged over in our effective description. However, there is an effect on 
the effective stress tensor measured at large distance from the brane, since the equation of 
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state (which is known only perturbatively near the GL threshold, or possibly numerically) 
is in general different than the one for smooth branes ()2.19|) . 

We can use these lumpy branes as the basis for the construction of lumpy blackfolds. 
It should be clear that their worldvolume is smooth, but the horizons of these blackfolds 
are inhomogeneous on the scale tq. The simplest example would be a lumpy black ring, 
built by bending a lumpy black string into a circular shape. This example serves also to 
illustrate a feature of lumpy blackfolds: the lumps in a rotating lumpy black ring will emit 
gravitational radiation, so the ring will lose mass and angular momentum as it evolves in 
time. In more generality, if the lumps on a blackfold extend along a direction in which 
the fluid velocity is non-zero, and if this direction is not an isometry, then the lumps 
moving along these orbits will give rise to a varying quadrupole and hence to gravitational 
radiation. 

Once this effect is taken into account, lumpy blackfolds are generically expected to 
evolve in time and not remain stationary. However, the time-scale for this evolution will 
be very long. The effect is only visible when the small scale is resolved, so it is suppressed 
by a power of ro/R. It will be further suppressed by the fact that gravitational radiation 
couples to a higher-multipole (quadrupole) and therefore is rather inefficient. In [29] this 
time scale was estimated for five-dimensional rings, and extending this estimate to black 
rings in D = 4 + n dimensions we find the time to be of order T gw ~ R(R/ro) n , longer by 
a factor (Z2/ro) n+1 than the short time-scale r^. Lumpy black branes may also be affected 
by GL-like instabilities, but these have not been investigated yet. 

3.5 Generalization to non-vacuum theories 

In (|2,ip and (|2.2p we assumed that the full dynamics at all wavelengths is described by 
vacuum General Relativity, i.e., the Einstein-Hilbert action with no matter nor cosmolog- 
ical constant. However, this is not actually necessary for our derivation of the equations of 
motion and the effective stress tensor. The only part of the field that is actually required 
to be governed by the Einstein-Hilbert action is the sector of short-wavelength degrees 
of freedom that we integrate out in order to obtain the effective stress tensor (|2.13|) . For 
the long-wavelength components we only require diffeomorphism invariance, which implies 
the equations of motion (|3,10p . Thus, the blackfold equations (|3.19p are enough to de- 
scribe neutral blackfolds in any configuration that, at small distances, is dominated by the 
Einstein-Hilbert term. For instance, this will be the case for blackfolds in the presence of 
a cosmological constant as long as 

r « \K\- 1 ' 2 (3.25) 

(see [13] for an explicit application), or for blackfolds in an external background gauge 
field as long as the typical length scale of the background field around the blackfold is 
much larger than tq. No restriction on R other than R^> ro needs to be imposed. 
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A different situation arises for charged blackfolds, since then the gauge field has short- 
wavelength components. The effective fluid is then charged, and additional current conser- 
vation equations must be added. This extension of our analysis will be discussed elsewhere. 



4 Stationary blackfolds 

Equilibrium configurations that remain stationary in time are of particular interest. For a 
blackfold, they correspond to stationary black holes. In this case it is possible to solve the 
blackfold equations explicitly for the worldvolume variables, namely the thickness tq and 
velocity u, so one is left only with the extrinsic equations for the worldvolume embedding 

We employ a general result proven in [30] for stationary fluid configurations: if dissi- 
pative effects must be absent, then the fluid (intrinsic) equations require that the velocity 
field be proportional to a worldvolume Killing field k = k a d a . That is, 

u = k/|k| (4.1) 

where 

|k| = V~labk a k b (4.2) 



and k satisfies the worldvolume Killing equatior 

D {a k b) = 0. (4.3) 

This does not necessarily mean that k is a Killing field of the background away from 
Wp+i, as the worldvolume could be at a locus of enhanced symmetry. However, this is not 
generic, and indeed when the blackfold thickness is small but non-zero, k should satisfy 
the Killing equations on a finite region around the worldvolume. Thus we assume the 
existence of a timelike Killing vector k^d^ in the background, 

V (/ A)=0, (4.4) 

whose pull-back onto the worldvolume determines the velocity field as in (|4.ip . The 
existence of a timelike Killing vector field is in fact a necessary assumption if we intend 
to describe stationary black holes. 

The contraction of the Killing equation (|4.4p with k^k v implies fc^9^|k| = 0, and using 
this it follows easily thaj^l 

^ = ^ln|k|. (4.5) 



For a stationary conformal fluid this is relaxed to the conformal Killing equation D^ a k b ) — \y a b- 
Conformal fluids do not dissipate the expansion of it since the bulk viscosity vanishes. 

n in [30] it is shown that the worldvolume fluid equations directly imply u b Dtu a — 9 a ln|k|, which 
implies the worldvolume projection of (|4.5[) . but not the orthogonal component of it, which will be used 
later below in eq. (|7.ip . 
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Since the expansion of u vanishes, the intrinsic blackfold equation (|3.22|) becomes 



<9 a ln|k| = d a lnr (4.6) 

so 

= constant . (4-7) 

In order to fix the proportionality constant in (|4.7p , we turn to the fact that the stationary 
blackfold describes a black brane with a Killing horizon. 

Above we have defined k as a Killing vector in the background. Its norm k 2 = — |k| 2 
computed with the worldvolume metric 7 a & is negative, i.e., it is a timelike vector. But 
recall that we regard the background metric as only one part, the far-zone metric at 
r 3> ro, of the full geometry. For r <C R the geometry is instead well approximated by 
the near-zone metric (|2.8|) . which matches with the far- zone metric in ro <C r -C R. It is 
natural to extend k as a Killing vector to all the geometry, including the near-zone (|2.8p . 
Using the metric gfu/° in this region the norm of k is 

5 (sho rt )^ r = ^ + r l UaUb \ k a k b = - (l - 3[\ |k| 2 (4.8) 

where we used (|4,ip and (|4.2p . Thus k becomes null as we approach the horizon r — > ro 
and indeed it is the null Killing generator of the horizon. Its surface gravity is easily 
obtained as 

^|k| 

k = -Li . 4.9 

2r 

Equation (|4.7p tells us that k is a constant over the worldvolume of the blackfold. That is, 
the surface gravity is uniform not only over the (n + l)-sphere but over the entire horizon. 

Eqs. (|4.ip and (|4.9p provide the general solution to the intrinsic equations for stationary 
blackfolds. To give them a more explicit expression it is convenient to choose a preferred set 
of orthogonal commuting vectors of the background geometry, ^ timelike and Xi spacelike, 
such that the Killing vector k is a linear combination of them 

k = £ + Er^' (4.io) 

i 

with constant f2j. Clearly i runs at most up to p. A convenient (but not necessary) 
choice is to take for £ the generator of asymptotic time translations, and for Xi the Cartan 
generators of asymptotic rotations with closed orbits of periodicity 2ir. They will often be 
Killing vectors themselves, but in principle only the linear combination in k need be so. 
Introduce a set of worldvolume functions R a (c) via their norms on the worldvolume, 



^0 = ^ 



IV, 



Ri = \x 



p+1 



w, 



(4.11) 



The R a must be regarded as part of the embedding coordinates X^(a). If £ is the 
canonically-normalized generator of asymptotic time translations then Rq is a redshift 
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factor between infinity and the blackfold worldvolume. The Ri are the proper radii of the 
orbits generated by x% along the worldvolume. The Oj are the horizon angular velocities 
relative to observers that follow orbits of £. 
The vectors 

Jhi e ' h = k Xi (412) 

(no sum in i) are orthonormal with respect to the metric j a b on the worldvolume. It will 
also be useful to regard them as vectors that extend into all of (|2.8p . 
Let us introduce the worldvolume spatial velocity field 

Vi{a) = ^% = ^4 (4-13) 



so that 



and 



dt R (a) 



1/2 



where 



RoVl ~ V 2 , (4.15) 



^ = £^ = ■4 ( 4 - 16 ) 



Thus |k| can be regarded as the relativistic Lorentz factor at a point in VV p +i, with a 
possible local redshift, all relative to the reference frame of ^-static observers. Plugging 
(|4.15p into (14. 9h we obtain that for given k and f2j the thickness ro is solved in terms of 
the R a as 

r Q {o) = ^^VHa). (4.17) 

The conditions that k and must remain uniform over the blackfold worldvolume 
were referred to in pQ as the blackness conditions and were imposed, invoking general 
theorems for stationary black holes — zeroth law of black hole mechanics and horizon 
rigidity [3U E21 E3] — , as necessary conditions for the regularity of the black hole horizon. 
Here instead we have derived them as general consequences of stationary fluid dynamics, 
where k and fij appear as integration constants. 

As a matter of fact it is also possible to work entirely within the framework of the 
effective theory and avoid any reference to the short-wavelength geometry of the horizon. 
Using only the fluid and thermodynamics equations one can derive, like in [30], that the 
variation of the local temperature T (|2.2ip along the worldvolume is dictated by the local 
redshift 

T = T (4.18) 

RoVT^V^ 
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The integration constant T can then be interpreted, using the thermodynamic first law 
that we derive below, as the overall temperature of the black hole. As expected, T = k/2tt. 
However, we think it is instructive to see how the concepts of black hole physics, like null 
horizon generators and surface gravity, are recovered in this scheme. 



5 Blackfolds with boundaries 

Typically branes (such as D-branes or cosmic strings) acquire boundaries when they in- 
tersect or end on other branes. This effect is accounted for by adding boundary terms to 
the effective stress tensor. However, black branes (and other fluid branes) may also have 
'free' boundaries without any boundary stresses. 

Consider the case where the worldvolume of the brane has a timelike boundary dW p +i. 
Assume there is a smooth (but otherwise arbitrary) extension W p +i of W p +i across the 
boundary. Introduce a level-set function f(cr a ) such that / > in W p +i and / < in 
Wp + i — Wp+i, with f(a a ) = on the boundary dWp+i. Then — d a f is a one-form normal 
to dWp+i pointing away from the fluid. The stress tensor is 

T ab = [{e + P)u a u b + P lab ]Q{f), (5.1) 

where is the step function. 

If the fluid is to remain within its bounds then the boundary must be advected with 
the fluid, i.e., it must be Lie-dragged by u, 

£u\ dWp+ J = (5.2) 

or equivalently, the velocity must remain parallel to the boundary, 

At the boundary, the stress-energy conservation equation (|3.13|) becomes 

((e + P)u a u b + P lab ) d a f\ QWp+i = . (5.4) 
Imposing (|5.3|) we find that the pressure must approach zero at the boundary, 

P \ow P+1 =°- ( 5 - 5 ) 

This is simply the Young-Laplace equation for a bounded fluid when there is no surface 
tension that could balance the fluid pressure at the boundary. That is, we do not introduce 
any such boundary stresses, although, as we mentioned they are of interest when one 
studies brane intersections. 

For a neutral blackfold, vanishing pressure at the boundary means 

r °!aw P+1 =0 < ( 5 - 6 ) 
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which has a nice geometric interpretation: the thickness of the horizon must approach 
zero size at the boundary, so the horizon closes off at the edge of the blackfold. 

If the blackfold is stationary, the condition (|5.6p means that |k| —* so the fluid 
velocity becomes null at the boundary. This may happen either because the boundary 
is an infmite-redshift surface, Ro — > 0, or perhaps more commonly, because the fluid 
approaches the speed of light at the boundary, 

An example of this was presented in [T], where a rigidly-rotating blackfold disk was con- 
structed and shown to reproduce accurately the properties of MP black holes with one 
ultraspin. 

The extrinsic equations (|3.20p must also be satisfied and they impose further con- 
straints on the worldvolume variables. In the stationary case, using (|4.5p and (|4.7p we see 
that if the extrinsic curvature of the worldvolume remains finite at dW p +\, then not only 
ro — > but also _L PM d^ro must vanish there at least as quickly as tq. 



6 Horizon geometry, mass and angular momenta 

The blackfold construction puts, on any point in the spatial section B p of W p +i, a (small) 
transverse sphere s n+1 with Schwarzschild radius ro(c). Thus the blackfold represents a 
black hole with a horizon geometry that is a product of B p and s n+1 — the product is 
warped since the radius of the s n+1 varies along B p . The null generators of the horizon 
are proportional to the velocity field u. 

If ro is non-zero everywhere on B p then the s n+1 are trivially fibered on B p and the 
horizon topology is 

(topology of B p ) x s n+1 . (6.1) 

However, we have seen that if B p has boundaries then ro will shrink to zero size at them, 
resulting in a non-trivial fibration and different topology. A simple but very relevant 
instance of this happens when B p is a topological p-ball. Then the horizon topology can 
easily be seen to be SP +n+1 = S D ~ 2 . 

To analyze the horizon geometry we go to the metric (|2.4p that locally describes the 
geometry of the blackfold to lowest order in ro/.R in the region r <C R- There we can 
choose a local orthonormal frame (dt,d z i) on the worldvolume, such that dt coincides in 
the overlap zone ro <C r <C R with the timelike unit normal to B p , 

= (d t f . (6.2) 

To lowest order in tq/R the worldvolume metric is Minkowski and the spatial metric on 
the horizon at r = ro is 

ds 2 H = (5ij + UiUj) dz l dz j + r^dtiL^ (6.3) 
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with itj = u ■ d z i. Then the local area density of the horizon an at a given point in B„ is 



„n+l 



a H = ^( n+ i)rg + y/l+Jj^Uj . (6.4) 
For a stationary blackfold the choice for dt and c^i was made in (I4.12D . so 

n" = ^-e, (6.5) 



and 



so the area density is 



= ; (6.6) 



0, -~ n+1 



oh 



Z (n+l) r O 



= n (n + l)(^Y i^ +1 K)(l-^(0) f • (6-7) 
The total area of the horizon is then 

An= [ dV (p) a u (a a ) , (6.8) 
J Bp 

where dV(p) denotes the volume form in B p . 

Again, we could have avoided any reference to the geometry of the horizon and worked 
instead exclusively with quantities defined in the effective fluid theory. The entropy density 
of the blackfold fluid is given in (|2.2U|) and after taking the relativistic Lorentz factor 
y/l — V 2 into account, the total entropy is 

s - L - (£ L «** « (i - w»t - % . 

(6.9) 

in agreement with the geometric area computed from (|6.8p . The geometric interpretation 
involves short-wavelength physics, but is useful for exhibiting how the blackfold construc- 
tion gives precise information about the entire horizon geometry, including the size of the 
s n+1 . 

The mass and angular momenta are conjugate to the generators of asymptotic time 
translations and rotations, which we assume are the vectors £ and Xi that we introduced 
in sec. I4"l 12 l Then 



M 



[ dV {p) T^n"? , Jt = - ( dV (p) T^ X V i ■ (6.10) 

J Br, JBr, 



Plugging here (|6.5p and the results from sec. 0] we obtain 

M = ^ (£)" L ^ " ^ (»+ 1 " V) , (Ml) 



12 In a more general case the relation may involve linear combinations, but this, although straightforward, 
comes at the expense of more cumbersome expressions. 
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and 

■* = T££(£)"" n ' L, dV <» ^-O^ ■ (6.12) 

It is easy to extract some interesting consequences of these results. Let us assume 
that all length scales along B p are ~ R and that the velocities and redshift are moderate 
(i.e., 1 — V 2 and i?o of order one) over almost all the blackfold. Then the two black hole 
length scales introduced in (jl.ip are 

£ M ~(r^ p )^3 £j~R, (6.13) 
and the small expansion parameter for the effective theory is 

(6.14) 



x D-3 

m \ fr_o_ 
~ \R 

It is interesting to compare the areas (i.e., entropies) of different blackfolds in a given 
dimension D. In [6] we introduced the dimensionless angular momentum j and dimen- 
sionless area a for a given masi^l 

a ~ 775^2- ( 6 - 15 ) 



M(GM)V(fl-3) £ M ' ,/>-2- 



The blackfold approximation requires j>l. Since ~ Tq i? p we find that 

a (j) ^j-s4=F . (6.16) 



This shows that for a given number of large non-zero angular momenta, the blackfold 
with smallest p is entropically preferred at fixed mass. This is just like we observed for 
single-spin MP black holes and black rings in [6]: for a given mass, the smaller p, the 
thicker the horizon, thus the cooler the black hole, and (since kAh ~ GM i.e., constant 
for fixed mass) the higher its entropy. 



7 Action principle and first law of stationary blackfolds 
7.1 Action for the embedding geometry 

We have presented in sec. U] the general solution to the intrinsic equations for a stationary 
blackfold. Given a Killing vector field k, eqs. (j4. 1|) . (|4,15p and (|4.17p allow to eliminate 
the variables ro(cr) and u a (a) in terms of the embedding functions X fM (a). The remaining 
extrinsic equation (|3.20[) that determines the embedding can be written, using (|4.5p . as 

K p =± Pfl <9 M In |k| n . (7.1) 
13 One should not confuse the dimensionless total area for fixed mass a with the blackfold area density 
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Using the result (|A.40|) from the appendix, this equation can be obtained from the action 

I[X fi (a)]= f d p+1 a^\k\ n (7.2) 

by considering variations of in directions transverse to the worldvolume. 

We can write it in a form that is particularly practical for obtaining the blackfold 
equations in specific calculations. First observe that the asymptotic time, conjugate to 
the vector £, is related to proper time t on the worldvolume by a factor of i?n- If we take 
the interval for the (trivial) integration over asymptotic time to have finite length (3, then 



/ d p+1 a >/=7|k| n = [ dV (p) R \k\ n . (7.3) 
JW p+ i Jt3p 

Using now (|4.15p we find 

=/3 / dV {p) R^\l-V 2 )^ 

= (5 jf dV {p) Ro(a) (^R 2 (a) -^^R 2 ^ . (7.4) 

We emphasize again that R a are among the worldvolume field variables X fl (a). Of course 
these enter as well through dV^ 

7.2 First law 

Using eqs. (j6.8|) . (|6.1ip . (|6.12p . it is straightforward to check that the action (|7.4p is 



I = i , { M -T, a ' J >-£a A ») 



(7.5) 



This identity holds for any embedding, not necessarily a solution to the extrinsic equations. 
Thus, if we regard M, Ji and Ah as functionals of X^^a), and consider variations at fixed 
surface gravity and angular velocities, we have 

SI _ ^ 5M ^ 5J, k 5A H 
^ = ° * SX^-^ ni SxT-^GJx^ = - (7 - 6) 

i 

Hence, solutions of the blackfold equations satisfy the 'equilibrium state' version of the 
first law of black hole mechanics. Conversely, the blackfold equations for stationary con- 
figurations can be obtained as the requirement that the first law be satisfied. If we regard 
k and f2j as Lagrange multipliers we may also say that blackfolds extremize the horizon 
area for given mass and angular momenta. 



14 In 1 it was asserted that the blackfold equations can be derived from the action 
J w d p+1 o y J~^T ah 'yab- This is not true in general, but is correct in the stationary case since T ahr y a b oc 
ro n «ik|'\ 
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In the Euclidean quantum gravity approach to black hole thermodynamics it is natural 
to take (5 to be the period of Euclidean time, j3 = 1/T. Using kAh/8itG = TS and eq. (17,51) 
we see that / is equal, up to a factor, to the Gibbs free energy G, 

0- 1 I = G = M-^VL i J i -TS. (7.7) 

i 

Therefore (|7,2p can be identified as the effective action that approximates, in the blackfold 
regime tq/R <C 1, the gravitational Euclidean action of the black hole [34] . 

It is also possible to find action functionals for general, possibly time-dependent black- 
folds by adapting the action principles developed for perfect fluids [35]. However, the 
usefulness of these actions, which must deal with constraints such as u 2 = — 1, appears to 
be somewhat limited so we do not dwell on them. 



8 Gregory-Laflamme and correlated thermodynamic insta- 
bility in blackfolds 

The blackfold approach must capture the perturbative dynamics of a black hole when the 
perturbation wavelength A is long, 

A»r . (8.1) 

These perturbations can be either intrinsic variations in the thickness ro and local velocity 
u, or extrinsic variations in the worldvolume embedding geometry X. In general, these two 
kinds of perturbations are coupled. A detailed analysis of the perturbations of solutions 
to the blackfold equations and their stability will be presented elsewhere. Here we extract 
some simple but important consequences for perturbations with wavelength 

r <A<i?, (8.2) 

i.e., perturbations for which the worldvolume looks essentially flat, ~ 0. In this case 

it is easy to see that the intrinsic and extrinsic perturbations decouple. 

It is instructive to perform the analysis for a general perfect fluid (|2.15p . and then 
particularize to the neutral blackfold fluid (|2.19|) . For simplicity we consider a fluid initially 
at rest u a = (1,0. . . ), with uniform equilibrium energy density e and pressure P. The 
flat worldvolume metric is parametrized, in 'static gauge', using orthonormal coordinates 
X° = t, X % = z\ i = l,...p and the transverse coordinates X m are held at constant 
values. Introduce small perturbations 

dP 

S £ , SP = —5e, 5u a = (0,v l ), 5X m = £ m , (8.3) 
ae 

and work to linearized order in them. The perturbed stress tensor is 

dP 

T u = e + Ss , T tl = (e + P)v l , T n = P + — 5s , (8.4) 
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and the extrinsic curvature 

5K ab m = d a d b C . (8.5) 
The extrinsic equations (|3.12p then become 

(edf+Pdf)r = 0. (8.6) 
Thus transverse, elastic oscillations of the brane propagate with speed 

4 = -f • (8-7) 

The intrinsic equations (|3.13p are 

d t T a + diT u = , d t T u + djT ji = , (8.8) 
which can be combined into 



For (ED we find 



d?T tt - dijT ij = . 



de 

so longitudinal, sound-mode oscillations of the fluid propagate with speed 

4 = -. (8.11) 

These derivations of eqs. (|8.7p and (|8.1ip are hardly new: they are conventional ways to 
obtain the speeds of elastic and sound waves — in fact they have been obtained in [361 EZ! 
for brane dynamics. They have a remarkable consequence: a brane with a worldvolume 
fluid equation of state such that 

PdP . 
> 8.12 

e as 

has 

44 <° ( 8 - 13 ) 

and so is unstable to either longitudinal or transverse oscillations with wavelengths in 
the range (18. 2p . For instance this happens in the simple case P = we with constant w, 
where the interpretation is easy (we assume e > 0): positive tension is required for elastic 
stability, but positive pressure is needed to prevent that the fluid clumps under any density 
perturbation. 

Neutral blackfolds have 

c - = -* = -^TI- (8 ' 14) 

and therefore are generically unstable to longitudinal sound-mode oscillations and stable 
to elastic oscillations in the range of wavelengths (18. 2D . 
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This instability is not unexpected. Black branes suffer from the Gregory-Laflamme 
instability [28], which makes the horizon radius vary as 

5r ~ e nt+ik ^ . (8.15) 

Here Q is positive real and thus the frequency is imaginary. The threshold mode for the 
instability, with S7 = and k = V k^k 1 ^ 0, has 'small' wavelength A = 2ir/k ~ r$ and 
therefore cannot be seen in the blackfold approximation. But the GL instability extends 
to arbitrarily small k, i.e., arbitrarily long wavelengths, and when k is very small it should 
be captured by the blackfold dynamics. 

The sound-mode instability corresponds precisely to this long-wavelength part, f2, k — > 

0, of the GL instability. Observe that sound waves in a blackfold produce 6e ~ 5P ~ 5r$ 

1. e., variations in the horizon thickness. Eq. (18,14p tells us that these are unstable, of the 
form (|8.15p with dispersion relation 

n = —L=k. (8.16) 

A simple inspection of the slope at the origin in figure 1 of [28] shows good numerical 
agreement with ()8.16j) . We leave a more precise derivation of this equation from a detailed 
GL-type analysis to future work. Note also that the dispersion relation (|8,16p indicates 
that the collective coordinate ro is a ghost (i.e., its effective Lagrangian — cj 2 {dt lnro) 2 + 
(d Zi lnro) 2 nas the 'wrong sign' for the kinetic term) F^l 

Moreover, observe that the Gibbs-Duhem relation dP = sdT, from (|2.17j) and (|2.18j) . 
implies in general that 

dP _ dT _ s 

de de c v 

where c v is the isovolumetric specific heat. Thus the black brane is dynamically unstable 
(to long-wavelength GL modes) if and only if it is locally thermodynamically unstable, 
c v < 0. This is precisely the content of the 'correlated stability conjecture' of Gubser and 
Mitra [38] . In fact our method gives a quantitative expression for the unstable dynamical 
frequency in terms of local thermodynamics as 



n = J-^-k, (8.18) 

V \ C v\ 

which as far as we know is a new result. 

The ordinary derivation of the GL instability involves a complicated analysis of lin- 
earized gravitational perturbations of a black brane and the numerical resolution of a 
boundary value problem for a differential equation (which is moreover compounded at 
small k since larger grids are required to avoid finite-size problems). Here we have shown 
that the long- wavelength component of the instability, (|8.16|) . can be obtained by an almost 



'The GL mode at threshold is instead tachyonic, since its dispersion relation has imaginary mass. 
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trivial calculation of the sound speed in a fluid I. In addition, the correlation between 
dynamical and thermodynamical stability follows as an elementary consequence of the 
thermodynamics of the effective fluid. In our opinion these results are striking evidence 
of the power of the blackfold approach. 



9 Closing remarks 

The formalism we have developed bears relation to two different earlier effective descrip- 
tions of black hole dynamics. The extrinsic part is a generalization to p-branes of the 
effective worldline formalism for small black holes, whose size tq is much smaller than the 
length scale R = I /(acceleration) of their trajectories or the wavelength of the gravita- 
tional radiation they emit |18|, [T9l [20] . The intrinsic part is similar to other fluid-dynamical 
formalisms for horizon fluctuations, such as the membrane paradigm and the fluid/ AdS- 
gravity correspondence. 

With respect to the first one, note that our formalism allows to consider time-dependent 
situations, which typically involve the emission of gravitational waves. This can be ob- 
tained by coupling the blackfold effective stress tensor to the quadrupole formula for 
gravitational radiation. This is also common in studies of gravitational wave emission 
from small black holes and from cosmic strings. However, accounting for the backreaction 
of this radiation on the blackfold requires going beyond the generalized-geodesic approx- 
imation and dealing with the notoriously subtle problem of the gravitational self-force 
[39]. 

To relate the fluid/ AdS-gravity correspondence to our approach take, instead of a 
neutral black brane, a near-extremal D3-brane. The blackfold formalism can be applied 
to it, toc0. The small scale corresponds to the charge-radius r q of the D3-brane (which 
is much larger than the non-extremality length ro), and in the overlap-zone r q <C r <C R 
where the blackfold effective stress tensor is computed, the metric is flat up to small 
corrections in r q /R. The blackfold method here could be regarded as an extension of the 
DBI approach to describe thermally excited worldvolumes. The difference with respect 
to neutral branes is that there is a region near the horizon where excitations with long 
wavelengths 3> r q are localized. One can take a limit, Maldacena's decoupling limit, 
to decouple all the far-zone effects from them. This region is asymptotic to AdSs x 
S 5 with radius r q , and far-zone effects are absent since they would give rise to non- 
normalizable modes in AdS and change the boundary geometry. Integrating the degrees of 
freedom in the asymptotically AdS region one gets only intrinsic, purely hydrodynamical 
collective modes. The effective stress tensor thus obtained is again of quasilocal type 
and in fact is the holographic stress tensor in AdS |40j . This is in principle different 

16 Observe that this is not a Jeans instability of the fluid (which has sometimes been suggested as possibly 
related to the GL instability) since gravitational forces within the fluid are entirely absent in our analysis. 
17 Charged blackfolds will be discussed elsewhere. 
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than the one in the blackfold approach (which is computed in an asymptotically flat 
overlap-zone) but clearly they are related. For instance, at the perfect fluid level they 
are the same. The main point is that in the fluid/AdS-gravity correspondence there is no 
extrinsic worldvolume dynamics, nor far-zone backreaction. But the less general dynamics 
comes with a number of advantages: first, it makes much simpler to compute higher- 
derivative corrections to the perfect-fluid dynamics. Second, the charge near extremality 
eliminates the sound-mode instability. Finally, this specific correspondence can be argued 
to describe the hydro dynamical regime of a strongly-coupled quantum Yang-Mills theory. 
For a generic blackfold it is not known whether there is a useful dual interpretation of that 
sort. 

Connections to the membrane paradigm are also suggestive but remain somewhat less 
precise. The membrane paradigm can be formulated for generic black holes, including 
vacuum ones, in terms of an effective fluid with a stress tensor of quasilocal type |41j . But 
in the membrane paradigm the boundary where the fluid lives is taken to lie right above the 
horizon, and not in an intermediate-asymptotic region (overlap zone) as in the blackfold 
approach. Perhaps the membrane-paradigm stress tensor and the one for blackfolds can 
be related through Komar integrals, at least in vacuum gravity. At any rate the membrane 
paradigm seems to capture only what we refer to as intrinsic fluid dynamics, and it is not 
clear to us how it could deal with the extrinsic embedding dynamics of a black brane. 
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A Geometry of embedded submanifolds 

We collect here several relevant definitions and results on the geometry of submanifold 
embeddings. Some aspects are more extensively discussed in [32]. 
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A.l Extrinsic curvature 

Assume the submanifold W is embedded as X^{a a ). The pull-back of the spacetime metric 
onto W is 

7ab = g^d a x»d b x» . (A.l) 
The first fundamental tensor of the surface is then 

h iw = ^daX^dbX" . (A.2) 

It follows easily that 

U\d a X v = d a X» , (A.3) 

and 

h\h\ = h% (A.4) 
so h^ v projects tensors onto directions tangent to W. Decomposing the metric as 

9^iu = h^ u -\- -L^ u , (A-5) 

we obtain the orthogonal projection tensor -L^ u , 

-Lp, d a X» = , = !_/ . (A.6) 

The extrinsic curvature tensor can be defined as 

V = h\h%Vxh" a = -h\h\V x ± p a ■ (A.7) 

The tangentiality of the first two indices and orthogonality of the last, 

±/K CTi / = L a v K vv ? = K p K ail v = (A.8) 

follows easily from this definition and the projector property (|A.4D , 

Following [T3], it is convenient to introduce the tangential covariant derivative 

V M = fc/V„ , (A.9) 

so 

K^ p = K a V„h a p . (A. 10) 

Applying the tangential derivative on (|A.4p one obtains 



Let v be any vector tangent to W. Then 



(A.ll) 



= ±V (A.12) 
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where 

v» = v u VvV IM . (A.13) 
Now let N be any vector orthogonal to W. Then 

N p K^P = NphSV^K" = -h u PV^N p . (A. 14) 

The symmetry 

K [H P = (A.15) 

follows from the integrability of the subspaces orthogonal to -L^ u . To prove this, assume 
that the latter is true, namely that there is a submanifold W defined by a set of equations 
/(j) (X) = such that df^ are a basis of one-forms normal to the submanifold. Any one- 
form normal to W is a linear combination of them so the subspace orthogonal to it is also 
integrable. It is always possible to choose a one-form N orthogonal to this subspace such 
that 

N p = drfiX) = V M /(X) , (A.16) 

so, using (|A.14|) . 

NpK^P = -K°h/X7 a X7 p f , (A.17) 

which is manifestly symmetric under fi *-* v. The converse statement that (IA.15|) implies 
the integrability of the orthogonal subspace, can also be proven by a straightforward 
application of Frobenius's theorem |42j . 

Background tensors t lil n 2 ... VlV2 '" can be pulled-back onto worldvolume tensors t aia2 ^ blb2 " 
using d a X^ as 

t aia2 .. blb2 ~ = d ai x^d a2 x^ ■ ■ ■ d h x ul d b2 x U2 ■ ■ ■ w.. 1 ' 1 " 2 - > ( A -!8) 

where 

8 b X u = 1 bc h up d c XP . (A. 19) 

Observe that even when t^ lP2 ,, l ' lV2 '" is a background tensor with indices parallel to W, 
in general V M i AtlAt2 ... lyiiy2 "' has both parallel and orthogonal components. The parallel pro- 
jection along all indices is related to the worldvolume covariant derivative D a t aia2 .. blb2 '" as 
in (|A.18p . This can be shown by using the equation that relates the connection coefficients 
for each metric, T pu and { a c fe } , 

d a X»d b X v h%T% = d c X° { c ah } - h%d a d b X p , (A.20) 

which can be proven by direct substitution of the definitions of each term involved. 
In particular, the divergences of tensors are related as 

#V---Vpi pw - = d ai X^ ■ ■ ■ D c t ca >- . (A.21) 

Such relations allow to dispense with the use of worldvolume coordinate tensors and deriva- 
tives in most formal manipulations. However, worldvolume coordinates are very practical 
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for explicit calculations and also allow us to highlight the distinction between intrinsic 
(parallel to W) and extrinsic (orthogonal to W) equations. 

Let us now consider the divergence of a totally antisymmetric tensor J (such as a 
current associated to a gauge form field) parallel to the worldvolume. It is easy to show 
that 

J-'„ 1 V / ,J W * 1 - = ( A - 22 ) 
holds as an identity. This implies that the conservation equation 

V M J wtl - = (A.23) 

is equivalent to the worldvolume conservation equation 

D a J aai -=0, (A.24) 

i.e., the orthogonal component of the current conservation equation (|A.23|) does not yield 
any additional equation. In particular, for a 1-form current one has 

= D a J a , (A.25) 

and continuity of charge is only meaningful as an intrinsic equation. This is in contrast 
to the conservation of the stress energy tensor, where the orthogonal component gives 
independent extrinsic equations (|3.12p . 

Let us now obtain more explicit expressions for the pull-back of the extrinsic curvature 
tensor onto W in terms of X )M (a), 

K ab P = d a X^ b X v K^P = -d a X»d b X v V^ u P . (A.26) 

The property (|A.6j) implies 

= d b X v V v (± a Pd a X°) = -K ab P + U p d b X»V u (d a X°) . (A.27) 

Expanding the covariant derivative in the last term and using d b X v d v = d b we find 

K ab p = v {d a d b x° + r^daXtdbX") , (a.28) 

which is reminiscent of the expression for the acceleration (deviation from self-parallel 
transport) of a curve — indeed (|A.12p makes this even more explicit. An alternative 
expression with this same feature can be obtained by performing some manipulations: 

K p d a d b x° = d a d b xp - hp a d a d b x a 

= d a d b xp - { a c b } d c xp + uW^rj, 

= D a d b XP + d a X»d b X u hP a Tl u . (A.29) 
In the second line we have used (|A.20p . Inserting the last expression into ()A.28P we find 

K ab p = D a 8 b XP + FP^daX^dbX" . (A.30) 
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A. 2 Variational calculus 

Consider a congruence of curves with tangent vector N, that intersect W orthogonally 

N» V = . ^ -V= ^ > ( A - 31 ) 

and Lie-drag W along these curves. The congruence is arbitrary, other than requiring it to 
be smooth in a finite neighbourhood of W, so this realizes arbitrary smooth deformations 
of the worldvolume X p -> + 

Consider now the Lie derivative of /i^ along N . In general, 

V = N^Vph^ + h pu VpN p + h w V v N p . (A.32) 

Then 

h lx x h v a £ N h Xa = hp X h u °N p V p h Xa + VV M iV" + hppVvNP . (A.33) 
The first term in the rhs is zero: 

h/h^N^Vphxa = -h^h^NPVp i_ A(T 

= hp x ±xa NPVpK* = , (A.34) 

and the other two terms can be rewritten using (|A. 14|) . so 

W = -\K Xh ^ £ Nhxa . (A.35) 

This implies 

N p KP = -\h^£ N h^ = ±=£ Nx /\h\, (A.36) 

2 v\ h \ 

where h = det hu, v ■ These equations generalize well-known expressions for the extrinsic 
curvature of a codimension-1 surface. The last one allows to derive the equation for a 
minimal- volume submanifold: 

Vol= [ yj\h\ => 5 N Vol = -y/\h\N p K p (A.37) 
Jw 

i.e., for variations along an arbitrary orthogonal direction N, minimal (actually, extremal) 
volume requires K p = 0. This is of course the variational principle for Nambu-Goto-Dirac 
branes. 

Consider now a more general functional 



/= / (A.38) 
Jw 

where $ is a worldvolume function. Then 

5 N I = £ N (v/W$) = s/\h\{-N p K p $ + N p d p $) . (A.39) 
Since N is an arbitrary orthogonal vector we have 

5 N I = & K<> =± pfi dftln® . (A.40) 
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